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Small organisms (e.g., bacteria) and artificial microswimmers move due to a combination of active 
swimming and passive Brownian motion. Considering a simplified linear three-sphere swimmer, we 
study how the swimmer size regulates the interplay between self-driven and diffusive behavior at 
low Reynolds number. Starting from the Kirkwood-Smoluchowski equation and its corresponding 
Langevin equation, we derive formulas for the orientation correlation time, the mean velocity and 
the mean square displacement in three space dimensions. The validity of the analytical results is 
illustrated through numerical simulations. Tuning the swimmer parameters to values that are typical 
of bacteria, we find three characteristic regimes: (i) Brownian motion at small times, (ii) quasi- 
ballistic behavior at intermediate time scales, and (iii) quasi-diffusive behavior at large times due to 
noise-induced rotation. Our analytical results can be useful for a better quantitative understanding 
of optimal foraging strategies in bacterial systems, and they can help to construct more efficient 
artificial microswimmers in fluctuating fluids. 
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I. INTRODUCTION 

Biological fl|, 0, Q and artificial microswimmers [l|, 0, 
H, S III] niove through a fluid by performing a series of self- 
induced shape changes [1, Q . Handicapped by their tiny 
size (typically a few micrometers for a bacterium [lo|). 
they are forced to swim at very low Reynolds numbers 
3? -C 1 [HI, [H, nil- Hence, in order to account for the 
resulting lack of inertia, the swimming strategies of mi- 
croorganisms are very different from those operative at 
human length scales. More precisely, since the fluid flow 
is reversible at low Reynolds number, locomotion in this 
regime is possible only if the swimming stroke violates 
certain time- reversal symmetries [1, [H, [IE ■ 

Stimulated by experimental advances [l|, 0, 0, S 
[Tsj , in recent years considerable progress was achieved in 
understanding the dynamics of deterministic microswim- 
mers models [H [HI, [H [H, [H, H [H, [13, H^. 
Yet, comparatively little is known quantitatively about 
the complex interplay between active self-motion, hy- 
drodynamic interactions, and thermal fluctuations in 
the surrounding fluid [H, [s^, [3l|, HI]. Very recently, 
first steps towards clarifying these issues were made by 
Howse et al. [6], who measured in their experiments the 
mean square displacement of chemically driven colloidal 
spheres, and by Lobaskin et al. [s^, who studied the 
Brownian motion of a triangular microswimmer at inter- 
mediate Reynolds numbers H \ hy combining Lattice 
Boltzmann simulations with a Langevin description of 
the swimmer in phase space. In the present paper, we 
would like to complement these investigations by concen- 
trating on the diffusive properties of mechanically driven 
microswimmers at low Reynolds numbers 3? <C 1. This 
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limit case is most relevant for bacterial motions and al- 
lows one to treat diffusive effects within configuration 
space. 

Specifically, we will focus on the following questions: 
How does the size of the swimmer affect its effective 
mobility in a noisy fluid? Which details govern the 
transition from quasi-ballistic self-motion to the diffusive 
regime? To shed light on these issues, we shall consider 
simplified quasi- linear p-sphere swimmers similar to those 
proposed by Najafi and Golestanian [iO]. More precisely, 
we will assume that internal forces, which generate the 
swimming strokes, are mediated by interaction poten- 
tials. This approach permits us to treat thermal diffu- 
sion effects within the Kirkwood-Smoluchowski scheme, 
originally developed to describe the diffusion of polymers 
in a fluctuating medium [sl, [s^, [S^] . 

Starting from the Kirkwood-Smoluchowski equation 
(KSE) ensures that hydrodynamic and stochastic forces 
are consistently coupled on the level of the Fokker-Planck 
description in configuration space [33|. Moreover, as dis- 
cussed in Sec. IIIH the corresponding Langevin equation 
can be used to derive closed analytical formulas for the 
orientation correlation time, the mean velocity and the 
mean square displacement of a single swimmer in three 
space dimensions. Although the analytical and numer- 
ical results in this paper refer to the case of a quasi- 
linear 3-sphere swimmer {p — 3), the formalism can be 
easily generalized to more complex models (e.g., fiexible 
p-sphere swimmer chains). Therefore, this approach can 
be generally very useful for studying Brownian motion ef- 
fects and hydrodynamic interactions in active biological 
systems at low Reynolds numbers. Furthermore, since 
it is straightforward to implement an external cofine- 
ment (e.g., tweezer or lattice potentials), the Kirkwood- 
Smoluchowski scheme can help to construct and optimize 
arrays [38l | of, e.g., micro-pumps that work efficiently on 
those scales where thermal fiuctuations in fluid become 
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non-neglible. 

Thus, purpose and content of the present paper can be 
summarized as follows: First, we will discuss a convenient 
formalism that allows to simulate active microswimming 
by means of Langevin equations and interaction poten- 
tials (Sec. [ni- Subsequently, we derive analytic results 
for the diffusion of a single swimmer fSec. thereby 
extending recent work of Golestanian and Ajdari [25| 
on deterministic swimmers. A thorough analysis of the 
single-swimmer case is instructive for a number of rea- 
sons: (i) Exact analytical results provide a useful test of 
numerical simulations, cf. Sec. lIVI (ii) Recently, Leoni et 
al. Q were able to experimentally realize an individual 3- 
sphere system similar to those considered here, (iii) Un- 
derstanding the noise-induced behavior of a single swim- 
mer is a prerequisite for understanding complex behav- 
ior and pattern formation in, e.g., self-assembling bacte- 
rial systems [l0|. (iv) Dynamical calculations as those 
presented below may provide a "microscopic" justifica- 
tion for purely probabilistic models of bacterial motil- 
ity t3^ - (v) Depending on the swimmer size, we find a 
rather sharp transition from purely Brownian to quasi- 
ballistic motions. From a (bio-)physical perspective, it is 
remarkable that the transition occurs when the 3-sphere 
swimmer model is tuned to bacterial parameters. Hence, 
loosely speaking, exploiting the interplay between Brow- 
nian motion and active swimming may indeed represent 
a useful strategy (not only) in nature. 



II. GENERAL THEORETICAL BACKGROUND 

We consider an ensemble of N micros wimmers, each 
consisting of p spheres. Neglecting inertia, the state 
of the system at time t is described by a set of coor- 
dinates (Xa) — {X(Qj)(i)}, where a = l,...,pN is a 
sphere index, and i — 1,2,3 labels the space dimension. 
Our subsequent analysis rests on the assumption that 
the stochastic dynamics of the swimmers in the fluid can 
be described, at least approximately, by the Kirkwood- 
Smoluchowski equation (KSE) [13, [11 HE HI, repre- 
senting evolution for the iV-particle probability density 
f{t, {x(^ai)})- We begin by recalling how the KSE can be 
translated into a Langevin equation for numerical simu- 
lations [IE 113] . Details of the swimming mechanism will 
be discussed in Sec. IIIBI 



A. Kirkwood-Smoluchowski and Langevin equation 

Considering a fluid of viscosity fj, and temperature T, 
the KSE reads 0] 

dtf = d^a,)^i^ai){f3j) { [d{l3j)U] f + ksJ d^pj)f} . (1) 

Here, ks denotes the Boltzmann constant, := 
d/dx(^ai): E^nd a sum is performed over equal double in- 
dices {ai). The (effective [13|) potential U governs all the 



internal and external swimmer interactions (see examples 
in Sec. Ill B")) . apart from the hydrodynamic interactions 
mediated by the fluid. The latter are included in the ten- 
sor Considering spherical particles of radius Oq, the 
'diagonal' components of J{ are given by 

5' 

^(m)(Qj) = — 7 7a = STT^Oa, (2a) 

where 6ij denotes the Kronecker symbol, and 7^ is the 
Stokes friction coefficient. The hydrodynamic interac- 
tions between different spheres are encoded in the 'off- 
diagonal' components ^(ai)(0j)j a ^ p. If these hydro- 
dynamic interactions are neglected, corresponding to the 
(infinitely dilute) limit case ^{ai){f3j) = 0, Eq. ([T|) reduces 
to an 'ordinary' Smoluchowski equation with a diffusion 
constant = ksT/ja for each sphere. 

Here, we are interested in the effects of hydrodynamic 
interactions, corresponding to ^(ai){i3j) 7^ 0. A sim- 
ple approximation for ^(ai){i3j)i obtained by solving the 
Stokes equation for a point-like source, is the Oseen- 
tensor [iolliH 

where ra/3i := Xai — xpi and Tap ■= \xa — ^c^j. How- 
ever, the associated diffusion tensor D'-' := /cbT JC*-' is 
not necessarily positive definite, leading to unphysical 
behavior if sphere separations become too small j42l . [i^ . 
In our numerical simulations we shall therefore use the 
improved approximation 

orM _arO , ("a + 4) ( x O rgfj^Tapj \ 

^(..)m - + 24^^ rip ~ r 

(2c) 

which was derived by Mazur [i^]. The additional term 
on the right-hand side of Eq. (|2c)) can be understood 
as the next-order correction in a radius-over-distance ex- 
pansion of the mobility tensor for two spheres [i^ . For 
spheres of equal size (oq = ap), the tensor IK = 'K^ de- 
fined by Eqs. ((2a|) and ll2cl) . reduces to the Rotne-Prager- 
Yamakawa tensor [H |46| |. While Eq. (Pc)) gives a more 
accurate description than Eq. (j2b|) at moderate densities, 
both expressions become invalid if the distance between 
spheres becomes very small. At very high densities, when 
sphere-sphere collisions dominate the dynamics, near- 
field hydrodynamics and lubrication effects must be mod- 
eled more carefully (47j . In the present paper, however, 
we focus on systems that can be described by Eq. (|^ . 

Unlike the Oseen tensor IK*^ from Eq. I|2b|). the tensor 
"K = IK*^ is positive definite for rap > aa + ap and thus 
can be Cholesky-decomposed in the form 

'^(ai)(P]) = lj^C(ai)(-ik)C{pj){^k)- (3) 

The decomposition ([3]) is crucial if one wishes to find 
a Langevin representation for the stochastic process 
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{^(aiil*)i_described by the KSE (P): Upon noting 
that la m 







(4) 



holds (for both Jf^ and M*^, one finds that the KSEjp 
corresponds to the foUowing Ito-Langevin equation [49| 

{kB7f'^C(^,)(^u)AB(-,k){t)- (5) 

Here, {^(/3j)} '■— {~'^(P3)U} comprises the deterministic 
forces acting on the spheres, and is a collec- 

tion of standard Wiener processes; i.e., the increments 
di3(^fc)(i) := i?^fc(t-|-di) — are independent Gaus- 
sian random numbers with distribution 

g-«V(2dt) 

J'{dS(^fc)W e [u,u + du]} = (-2^^^) 1/2 



(6) 



and, according to the Ito scheme [43], the coefficients 
C(ai){-ik) are to be evaluated at time t. For complete- 
ness, we still note that, upon formally dividing by At, 
the stochastic differential equation ([5]) can be rewritten 
as a 'standard' Langevin equation: 



(7a) 



where X^^^-^lt) := dX(ai){t)/dt denotes the velocity, and 
^{■yk){t) d_B(^j.) (t)/dt represents Gaussian white noise. 



I.e. 



(e(m)(i)) 



0, 

SapSij S{t - t'). 



(7b) 
(7c) 



B. Swimming mechanism 



Having discussed how to implement fluctuations, we 
still need to specify the swimming mechanism. To this 
end, consider two spheres a and P that form the leg of 
a swimmer. We assume that the internal forces between 
beads a and f3, which generate the swimming stroke, can 
be derived from a time-dependent potential of the form 

[/iog(i, dafj) = Y {dafs -[£ + \ap sin(cji -I- ipap)]Y, (8a) 

with c?ct/3(t) := \X a{t) — X fj{t)\ denoting the distance be- 
tween the spheres, Xa/i the approximate amplitude of the 
stroke, and t > aa+ap + Xap the mean length of the leg. 
The potential U\cg gives rise to two characteristic time- 
scales: the driving period :— 27r/cj and, for a sphere 
of mass Ma, the oscillator period Tq := 2tt/ y^^o/M^. 
Since we are interested in the over-damped regime, these 
time scales must be long compared to the characteristic 
damping time := Ma/ja. More precisely, we have to 
impose 
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(8b) 



corresponding to slow driving and fast relaxation. The 
constraint (|8b[) ensures that our swimmers behave similar 
to a shape-driven swimmer [2^ . 

In the remainder, we shall focus on 3-sphere swimmers, 
representing the smallest self-swimming system within 
our approach (two-sphere swimmers can achieve active 
locomotion only due to collective effects [l3|)- We con- 
sider three spheres (a, (3, 7) forming a swimmer with 
central sphere /3, e.g., in the case of a single swimmer 
(a,/3, 7) = (1,2,3) with middle sphere /3 = 2. The legs 
are given by d^p := Xp — X^ and dp^ := X^ — Xp, 
and we still define normalized connectors 



dap/dap , n 
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dpj/dp^ 



In order to ensure that the swimmer moves quasi- 
linearly p5l | , we introduce a stiffness potential 



lin ("a/3 • np^ - i f , 



(9) 



which, for K ^ fcoi penalizes bending. The resulting 
force components ^'(q") ■— — c^(ai)C^iin read explicitly 

Q 



^ak 



dap 

Q 



{S^k 



T,p%naf)k)np^i, (10a) 



^-ik — ~i — {Sik — npjinpjk) napi, 
dpj 



^I3k 



( plin I 7Tilin\ 



where Uapi := dapi/dap, and 

Q := Kf {uap ■ np^ - 1) . 



(10b) 
(10c) 

(lOd) 



Equations ([5]) and ^ provide a convenient way of mod- 
eling and simulating rigid or flexible p-sphere swimmers 
by means of potentials. We note that, by construction, 
the sum over the internal swimmer forces is zero. The 
total potential U appearing in the KSE ^ is obtained 
by summing over all effective interaction potentials ([8]) 
and (El). 



III. ANALYTICAL RESULTS 

We next summarize formulae for the correlation time of 
the orientation vector, the mean swimmer velocity and 
the spatial mean square displacement of an isolated 3- 
sphere swimmer. These analytical results can be ob- 
tained from the Langevin equation ([5]) by using the Oseen 
approximation !K ~ and their explicit derivation is 
discussed in the Appendix [XI 

A swimmer's motion can be characterized by its geo- 
metric center 
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R{t) := -(X1+X2 + X3 



and the orientation vector 



Nit) 



X3 — Xi 



IX, 



^ll 



(11a) 



(lib) 
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We are interested in determining the mean square dis- 
placement 

Vn{t):={[R{t)-R{0)]^), (12a) 
and the correlation function 



2)^(i) := {N{t)N{Q)), 



(12b) 



where the average is taken over fluctuations in the fluid 
(i.e., over all realizations of the Wiener process). 

It is convenient to discuss Dff{t) first. For a de- 
terministic initial state N{0) = {Nk{0)), we can write 
"^Nit) = {Nk{t)) Nk{0) with a summation over equal in- 
dices. To obtain an analytical formula for {Nk{t)), we 
assume that Eq. ((8bl) holds true and that bending is 
neglible, K ko. Then the swimmer behaves like a 
stiff, shape-driven Najafi-Golestanian [23| swimmer and 
we can approximate 



di2 ■■= 


X2 




(13a) 


d23 ■ = 


X3 


- X2 ~ AT ^23, 


(13b) 


di3 := 


X3 


-Xi ~ AT (di2 -f ^23), 


(13c) 


where, cf. Eq. ([Sa 


1, 






di2 


= e 


-1- A12 sin(tjt -t- ipi2), 


(13d) 


d23 


= e 


+ X23 sin(cji 4- 9323)- 


(13e) 



Adopting the Oseen approximation !K ~ , one can de- 
rive from the Langevin equations (O the following linear 
evolution equation (see App. lA 1[) 



(Nk) 



kR7 
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[1 


{di3_ 




- 1 




{, ai 


03 / 



(Nk). (14) 



where (Nk) (diV^ (i)/dt} . The l/a^-parts are contri- 
butions to the rotation rate due to noise on the spheres, 
whereas the 1 / 3-contribution is a correction due to hy- 



drodynamic interactions. Equation (I14|) can be solved 
exactly. The solution exhibits an exponentially decaying 
oscillatory behavior due to the periodicity of the swim- 
ming stroke di3. However, for £ ^ maxjAa/j} it suffices 
to approximate di3 ~ 2^, yielding an exponential decay 
TlNit) ~ exp(— t/rjv) with orientation correlation time 



TN 



[ kB7 
\ 167r^^3 



ai 



03 



1 



(15) 



The time parameter tn not only determines the tem- 
poral correlation of the orientation vector, it also plays 
an important role for the dynamics of the geometric cen- 
ter R{t). As shown in App. lA 2 a[ the mean swimmer 
velocity {R(t)) is governed by the equation 



where 



(^(ifc)) 

{Fi3k)) 



B3 {dl2k) + C {d23k) 



C{d 



12k 



- C2 
^ ' Bi {d23k) 



B1B3 - C2 



(16a) 

(16b) 
(16c) 
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Att^i \di3 
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— (- 
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2tt^ di2 ' 


B3 
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1 
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72 


73 


27TfI d23 ' 


C 


1 

72 


1 

47r/i 


(— + — 

\ dl2 d23 



are the noise-averaged internal forces on the first and 
third sphere (the force on the central sphere can be elim- 
inated by virtue of -F(ifc) +^"(2/0) +-F'(3fc) = 0), respectively, 
and 

{di2k) = (Nk) du + (Nk) di2, (16d) 
{d23k} = {Nk) d23 + {Nk) d23 (16e) 

the mean change of the leg vectors due stochastic rota- 
tions and swimming strokes, with abbreviations 

1 



(16g) 
(16h) 
(16i) 

(16j) 



Since the quantities {Nk), {Nk), dap, and dap are known, 
Eqs. (I16p provide a closed analytical result for the mean 
velocity {Rk) of a shape-driven swimmer (within the Os- 
een approximation). In particular, Eqs. ()16p general- 
ize the corresponding velocity formulas for a determin- 
istic swimmer, recently obtained by Golestanian and Aj- 
dari psj . to the "noisy swimming" regime. 

For realistic swimmer parameters the orientation cor- 
relation time Tjv is typically much larger than the driving 
period = 27r/tt'. In this case, the rather lengthy re- 
sult p6)) can be considerably simplified (see last part of 
App. lA 2"al for details) to read 



{R{t))=V{N{t)). 



(17a) 



Here, V denotes the stroke-averaged velocity (i.e, over 
an interval \t,t -l-T^Jf) of the corresponding deterministic 
swimmer |25| . By using the approximation (|17ap instead 
of the exact results one neglects mean velocity oscil- 
lations on small time scales. For example, when consider- 
ing equal-sized beads with = a and £ ^ max{a, Aq^}, 
then 



V 



7 

24' 



A12A: 



23 



sin A(/j, 



(17b) 



where At/? := ipi2 — <^23 is the phase difference of the 
leg contractions, and higher order terms in a and \ap 
are neglected. Integrating Eq. p?|) with {N{t)) ~ 
Ar(0) exp(— i/rjv), we obtain for the position mean value 
of the swimmer the simple approximate result 

{R{t)) ~ R{0) + Vrjv [1 - exp(-t/Tw)] N{0); (18) 

i.e., in the asymptotic limit t 00, 

{R{oo)) = R{Q) + V tnN{0). (19) 
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Finally, let us still consider the mean square displace- 
ment Dfl(t) := {[R{t) - -R(0)]2) for a stiff, shape-driven 
3-sphere swimmer described by Eqs. (fT3| . As discussed 
in App. lA 2 bl by starting from the Langevin equation 
for R{t), one can show that T)ji{t) decomposes into the 
form 



(20a) 



where the first part 



1 fceT 



9 TT/i 
9 TT/i 
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— h 
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ai 


0,2 


0.3 


rt 






ds 




+ - 


3 




( 



1 

^23 



Z^ + ^ + Z^ (20b) 



1 

di3 



comprises passive Brownian motion contributions due to 
thermal diffusion of the spheres (first line) and hydro- 
dynamic Oseen interactions between them (second line), 
while the second part 



ds 



du {N{s)N{u)) (20c) 



is the contribution due to active self-swimming. Similar 
to Eq. pT]) . the expression l|20cp is valid if the orienta- 
tion correlation time is much larger than the stroke pe- 
riod. Tat T^. Inserting the above result {N{t)N{s)) ~ 
exp(— t/r^r) and approximating di2 = d23 = di3/2 ~ i, 
we obtain for the spatial mean square displacement 



2)fl(i) 



1 fcnT 



1 1 1 5\ 

— + — + — + - ]t + 

oi 02 as I 



9 TT/i 

2V^TN [i + rAr(e^*/"" 



!)]• 



(21) 



The approximate result (j2ip provides a coarse-grained 
stroke-averaged description of the translational diffusion 
(details of the swimming stroke are encoded in V). An 
analogous formula can be used to describe the diffusion 
of a spherical, chemically driven microswimmer Q. By 
virtue of Eq. (|2ip . we can readily distinguish three dis- 
tinct regimes: 

(i) For t tn we can expand the exponential term to 
linear order and find 



(22a) 



i.e., passive Brownian motion dominates on very short 
time scales. 

(ii) For t < Tj^ we need to include terms quadratic in 
t/TN and obtain 



'Dnit) ^ D^(t) + 



TN 



(22b) 



i.e., ballistic motion can dominate on intermediate time 
scales provided V'^/tn is large enough (cf. examples be- 
low). 



(iii) For t^r^, we recover diffusive behavior 



lim 

t — >oo 



(22c) 



If the swimmer is constructed such that 2V'^tn ^ 
D^j^{t)/t, then its diffusive behavior at large times is due 
to noise-induced rotation (with persistence time r^r). 

In a sense, the above results provide a 'micro- 
scopic' justification for the assumptions made by Lovely 
and Dahlquist [39], who studied purely probabilistic 
models of bacterial motion. We also note that the asymp- 
totic behavior for f 3> tat and 2V^tn ^ D^{t)/t agrees 
qualitatively with results reported by Lobaskin et al. [s^l 
for triangular swimmers in the moderate Reynolds num- 
ber regime 3? ~ 1. In this context, we also mention recent 
work by Golestanian et al. [lOl , who discuss similar scal- 
ing relations for the diffusion of phoretic swimmers. 

In the remainder, we are going to compare the ana- 
lytical predictions with results of computer simulations, 
based on a direct numerical integration of the Langevin 
equations (O for the spheres. Our main focus is on the 
transition from the passive Brownian motions to the ac- 
tive swimming regime. 



IV. SWIMMER TUNING & NUMERICAL 
SIMULATIONS 

We simulate a single 3-sphere swimmer described by 
the interaction potentials ([8]) and ([9]) and governed by 
the Langevin equation ([5]) . We consider identical spheres 
of radius = a, mass Ma = M , and equal stroke ampli- 
tudes Ai2 = A23 = A. The density of the spheres is chosen 
as p = lO'^ 'kg/ve? (water), and the fluid is water at room 
temperature [n = 10"^ kg/(ms), T = 300 K]. By fixing 
the parameters of the spring and bending potentials as 
kfi = 10-4 kg/s^, uj = 10^ Hz and K = lOfco, we sat- 
isfy the time scale condition (jSbp while ensuring that the 
swimmer behaves approximately stiff and shape-driven. 
The velocity V of the swimmer is optimized by choosing 
= 7r/2, cf. Eq. P5a|) . 

We are primarily interested in understanding how a 
change of the swimmer size may affect the diffusive be- 
havior. We therefore fix the ratios i/a = l/X = IQ and 
only vary the leg length t from 1 to 10 /xm in our simu- 
lations (i.e., the mean swimmer length is A = 2t). Put 
differently, we scale the swimmer proportionally by vary- 
ing I. Having specified all parameters, it is useful to 
summarize the relevant formulae for our choice: 
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— f 

24 



sinAifi =0.292 



£ 



kal \ 3a 



= 0.473 



1 fcsT /3 
a 



9 TT/i 



= 5.127 



/im^ 



(23a) 
(23b) 
(23c) 
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FIG. 1; (Color online) Orientation correlation function for 
four different swimmer sizes 2£. Symbols represent aver- 
ages over 100 trajectories, numerically calculated from the 
Langevin equation ([SJ using parameters as described in the 
text. The dashed lines depict the theoretically predicted ex- 
ponential decay I'jv(t) — exp(— t/rjv) with orientation cor- 
relation time rjv determined by Eq. (llSp . It is remarkable 
that changing the swimmer size by one order of magnitude 
increases the correlation time by three orders of magnitude. 

and 

£5 

2VVAr = 0.080 5-. (23d) 

It is remarkable that increasing the swimmer size by one 
order of magnitude increases the orientation correlation 
time Tat by three orders of magnitude. 

Figures [1] and [2] depict the results of numerical sim- 
ulations (symbols) of the Langevin equations (O and 
also the corresponding theoretical predictions (dashed 
lines). The numerical data points represent averages 
over 100 trajectories with identical initial conditions. 
More precisely, at time t = Q the swimmer is point- 
ing along the xa-axis with the first sphere being located 
at the origin, i.e., Xi(0) = 0, ^2(0) = ^iV(O), and 
^3(0) = [2£- Asin(A(^)]Ar(0), where N{Q) = (0,0,1). 

As evident from the diagrams in Fig. [1] and [H the re- 
sults of the numerical simulations are very well matched 
by the theoretical curves over several orders of magni- 
tude in time. In particular, for a leg length in the range 
£ ~ 5 /im (red /green "x "-symbols) one readily ob- 
serves the three aforementioned regimes: (i) Brownian 
diffusion at small time scales t <^tm, (ii) ballistic behav- 
ior at intermediate time scales, and (iii) quasi-diffusive 
behavior due to noise- induced rotation for t 3> t^t. 

We conclude the discussion of the numerical results by 
addressing a few technical aspects that might be rele- 
vant and helpful for future simulations. When consid- 
ering ensembles with > 1 swimmers the computa- 
tionally most expensive step is the Cholesky decompo- 
sition of the diffusion tensor, see Eq. ([3]), which is ap- 
proximately of the order 0[{pNY] [3g|. It is also worth- 
while to briefly comment on the choice of the time step 



FIG. 2; (Color online) Mean square displacement T>ii(t) :— 
{[Rit) - R{0)f} divided by time t for the same set of parame- 
ters as in Fig[T] The dashed lines correspond to the analytical 
formula (|21|l . The dynamics of small swimmers {£ — 1 fim, 
blue circles) is dominated by Brownian motion on all time 
scales, whereas big swimmers {£ = 10 fim, black diamonds) 
can move ballistically for several minutes. In the interme- 
diate region (£ ~ 5 fim, red "-f" /green "x") we observe a 
ballistic transition from ordinary Brownian motion at small 
times t <^ tn to noise-induced rotational diffusion at large 
times t ^ tn- Since rjv ~ the transition from Brownian 
to quasi-ballistic motion is very sharp. Interesti ngly enough, 
typical sizes of bacteria lie in or near this niche [lOl ] . 

dt in the Langevin simulations. Ideally, one would like 
to choose dt smaller than the smallest dynamical time 
scale in the system, which for our model is given by the 
damping time — M/(6TTfia). For the swimmer param- 
eters considered here we find ~ 2 x 10^^ (£//im)^s, 
which means that adopting dt would not allow us 

to simulate experimentally accessible time scales in the 
seconds range. Since we are not interested in the dy- 
namical details at very short times scales, we choose in 
our simulations the time step larger than Tj, but much 
smaller than the period of a swimming stroke by fixing 
dt = lO^^T^ for £ < 5Aim and dt lO^^T^ if ^ > 5^ni. 
We verified, however, that for intermediate time scales 
(of the order of a few stroke periods 6 x 10~^ s) the 
numerical results for the mean square displacement and 
other statistical observables agree with those obtained for 
very small time steps dt = O.IT.^. Generally, a satisfat- 
cory resolution of the bending and relaxation dynamics 
of the legs/spheres would require dt < {M/Ky/^ and 
dt < (M/A:o)i/^, respectively. 

V. CONCLUSIONS 

Understanding the interplay between Brownian mo- 
tion, hydrodynamic interactions, and self-propulsion is a 
prerequisite for understanding the dynamics of bacteria 
and artificial swimming devices at the microscale. In the 
first part, we discussed how one can model these phe- 
nomena by means of stochastic processes (overdamped 
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Langevin equations). Subsequently, as a first appli- 
cation, we focussed on the size-dependence of diffu- 
sive behavior at low Reynolds numbers (31 ^ f) for 
a quasi-linear 3-sphere swimmer model. Our theoreti- 
cal analysis complements a recent experimental study by 
Howse et al. , who investigated the diffusion of chem- 
ically driven, spherical colloids, and theoretical work by 
Lobaskin et al. [ssj . who considered the Brownian dy- 
namics of an artificial triangular microswimmer at mod- 
erate Reynolds numbers (3? ~ 1). 

Starting^from the Kirkwood-Smoluchowski equa- 
tion [sl, m, m, [131, we derived analytical results for 
the orientation correlation time, the mean velocity, and 
the mean square displacement of an overdamped, quasi- 
linear 3-sphere swimmer [20l . [25| . Analytical formulae 
as derived here are useful for testing numerical simu- 
lations and (in) validating simplified probabilistic mod- 
els [1^. Moreover, they provide detailed insight into the 
size-regulated transition from predominantly random to 
quasi-ballistic motions. 

The proposed method of modeling swimmers by effec- 
tive potentials within a Langevin scheme can be readily 
extended to study complex behavior in larger swimmer 
ensembles. However, at high swimmer densities, colli- 
sions and near-field hydrodynamics affect the diffusive 
behavior [13, Ull, [s^l and it will be necessary to modify 
the hydrodynamic interaction tensor accordingly. Gen- 
erally, a useful dimensionless quantifier for the efficiency 
of active swimming relative to diffusion is given by the 
ratio 



A/V" 



(24) 



The denominator corresponds to the time needed by a 
unperturbed swimmer of velocity V to move one body 
length A, and the swimmer geometry is encoded in the 
orientation correlation time rjv. For e f (e <C f ) self- 
propulsion is effective (non-effective). For non-isolated 
swimmers, r^r is not only determined by thermal effects, 
but also by collisions with other swimmers [5l|, . 

With regard to future studies we note that the com- 
bination of thermal fluctuations and hydrodynamic cou- 
pling might also lead to interesting behavior in simple 
arrangements of microswimmers. For example, experi- 
ments have shown that colloidal spheres localized in an 
array of optical trap s create memory effects (s^l and 
driven vibrations [54| . Our formalism provides a starting 
point for the investigation of many self-propelled bod- 
ies in separate potentials, which could be helpful for in- 
terpreting experimental data of trapped bacteria or for 
constructing pumps from a collection of microswimmers. 

To summarize, the above results may provide guid- 
ance for constructing artificial microswimmers or pumps 
that work efficiently in the critical transition region that 
separates Brownian from quasi-deterministic motions. 
In particular, by tuning the parameters to the narrow 
cross-over region one could construct swimmers that ex- 
plore with high probability a maximized volume fraction 



within a given period of time. The fact that many bac- 
teria live near or exactly in this niche [l3| suggests that 
this may indeed be a useful strategy. Thus, exploiting the 
interplay between noise and active self-motion could lead 
to novel applications [55| . e.g., with regard to the con- 
trolled transport [s^ and distribution of chemical and 
biological substances in small scale technical devices or 
even within the human body. 
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APPENDIX A: CALCULATIONS 

This appendix provides derivations of analytical results 
for the orientation correlation function, the mean veloc- 
ity, and the spatial mean square displacement from the 
Langevin equation Our calculations are based on the 
following simplifying assumptions: 

(i) The motion of the 3-sphere swimmer is approxi- 
mately stiff and shape-driven, i.e.. 



where 



di2 ■■= 


X2- 


Xi 


-Ndi2, 


d23 := 


X3- 


X2 


~Nd23, 


di3 := 




Xi 


~ N {di2 + d23) 


di2 


= i + 


A12 




d23 


= 1 + 


A23 


sm[ijjt — A(p). 



(Ala) 
(Alb) 
(Ale) 



(Aid) 
(Ale) 

In this case the internal forces, which generate the swim- 
ming strokes, point along the swimmer's axis and we may 
write 



Fu =giN,, F3, = 53 N, 
and, with Fai + F2i + Fu = 0, 

F2^ = -{gi+g3)N„ 



(A2a) 



(A2b) 



where ga is the force amplitude, and Ni denotes a com- 
ponent of the orientation vector 



Nit) := 



Xi 



1x3 - Xi 



(A2c) 



(ii) We adopt the Oseen approximation, i.e., IK ~ 
where 



o 



1 



87r/i dap 
1 



5u + 



daPidajSj 



The second line follows from assumption (IA1[) . 



(A3) 
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1. Orientation correlation function 

Given a deterministic initial state Ai^(O) = {Nk{0)), we 
are interested in DAr(t) := (iV(t)7V(0)) = {Nk{t)} ^fe(O). 

The first step is to find the stochastic differential equa- 
tion (SDE) for the orientation vector N(t). This can be 
achieved by virtue of the Ito formula [43, l57| , yielding 

dNkit) = 9(„,)iVfcdX(„,) + 

D(ai){f3f)d(ai)d{l3])Nkdt, (A4) 

where dX/^ai) (0 is governed by Eq. ^ , and the diffusion 
tensor is given by D := fcsT IK ~ fc^T !K'^ . 

To determine {Nk{t)), we use that 
(C'(ai)(7fe) di?(^fc)(t)) = for an Ito SDE, and there- 
fore 

d(^(Qi)(t)) = {Ji{ai){i3j)F(j3j))'it. (A5) 

Taking the average of Eq. (jA4p and inserting (|A5p . we 
obtain 

{Nkit)) = (ai)(P])F(l3j)d(c,i)Nk) + 

{D(ai)(P])d{ai)d(i3j)Nk) , (A6) 

where d{Nk{t)) := {dNk{t)/dt). We next evaluate the 
two terms the rhs. of Eq. (IA6p separately. To this end 
we note that 

d(<.j)Nk = -^{hc - Sia)iSjk ~ NjNk), (A7) 

"13 

yielding for the first term 



m 

di3 



i6^k - N,Nk). (A8) 



Using Oseen approximation (jA3p . we find 

91 (33 + 51) 



W(3i)(/3j)-F'(/3i) 



53 

73 



N,.. 



and, similarly, oc Ni. Taking into account 

that 



N,{d,k-N,Nk) 



the first term on the rhs. of Eq 
Eq. (|A6p reduces to 



0, (A9) 
6]) vanishes, and 



{Nk{t)) = (L'(„,)(0j)a(„,)9(;3j)^fe) • 

By virtue of Eq. (|A7p , we find 
1 



(AlO) 



"13 



(53/3 - (5l/3)((53a - <^la) X 



To obtain the rhs. of (jAlOp . we still need to contract 
with the diffusion tensor Di^^i-^(^pjy This results in two 
contributions: From the diagonal part we get 



7a 



1 



{SNjNkNj - Nj6jk - NkSjj - Nj6jk) 



1 ^ 1 \ 2iVfc 
.71 73/ rf?3 ' 
while the hydrodynamic off-diagonal terms give 



(Al2a) 



ksT ^ {I - Sal3){S3a - Sia)iS3i3 ~ Sip) : 
+ N,N,) 

1 

87r/i da/? 
"13 



Combining the two contributions we find 



{Mt)) 









-kB7 (— - 








13 J 


d\3 



2nfi dl 



13 



(A12b) 



(Nk), 



{3NjNkN, - N,5,k - NkS,, - N,5,k) ■ (All) 



which for 7q, — 6niiaa gives Eq. (jl4p . 

2. Motion of the geometric center 

Following a similar procedure, we can derive analytical 
expression for the mean velocity {R(t)) and the mean 
square displacement 'Dji{t). 



a. Mean velocity 

First, we would like to determine the mean velocity 
{R{t)) of the swimmer's geometric center 

R{t) i(Xi + X2 + X3). (A13) 

Averaging the stochastic differential equation 

dR{t) := i(dXi+dX2+dX3) (A14) 

o 

with respect to the underlying Wiener process and divid- 
ing by dt, we obtain 

1 ^ 

(Riit)) = o ) 

a = l 

= ^E(%^ + (1-M5^(«)(«)%.))- (A15) 

a = l \ / 
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Considering as before a stiff, shape-driven swimmer and 
Oseen interactions = IK"^, we have i^(/3i) — ff/j-ZV^ and, 
therefore. 



(1 - Safi)^{a^){|3j)F(Pj) 



E 



^ AtT^J. dap ' 



(A16) 



Inserting this into Eq. (jAlSp gives 

1 
3 



7i 72 



.92 



47r/x di2 
91 

51 



73 
53 



47r^ di3 

93 
A-KfJ. d23 

92 



^TTfi dl3 A-K^JL ^23 , 

Since the internal swimming forces sum to zero, we may 
ehminate g2 by using g2 = —{gi + .93) to obtain 



{R^{t)) 



{N^gl) 
{N^gz) 



1 

71 
1 

73 



1 

72 
1 

72 







--)] 


47r^ 




d23/. 


±- 








Kdi3 





4^ (N.gi) + ^ (A^.53) 



(A17) 



Consequently, to find {Ri{t)), we still need to deter- 
mine the mean forces = {Nigp) on the first and 
last sphere, (3 — 1,3. This can be achieved as follows: 
Equation (|Aip implies that 



(di2fc) = {Nk)di2 + {Nk)di2. (Al8a) 

fe) = {Nk)d23 + {Nk)d23- (A18b) 

On the other hand, from the definition of the vectors dap 
and the Langcvin equations for dXa, we have 

(rfl2fe) = {'K(2k)(P])F{pj) ~ 'Ki^-^k)(l3j)F(i3j)) 
{d23k) = (J^, 



(Al9a) 



'•(3fc)(/3j)^(/3j) 



^{2k)if3j)F(^f)j})- (A19b) 



Inserting the explicit expressions for Ji(ai)(i3j) and F^pj-^, 
Eqs. (jA19|) can be rewritten as 



{d- 



12k 



1 

71 



1 



1 



72 27r/x di2 

1 1 1 1 

72 47r^ \di2 d23 
-Si (A^fc5i) - C {Nkg3) 



(Nkgi) - 



1 

di3 



and 

{d23k) 



1 

72 

1 

72 



1 



47r/i 
1 

73 27r^ d23 



1 

rfl2 
1 



1 

^23 



1 

dl3 



{Nkg3) 



C{Nkgi)+B3{Nkg3). 



{Nkg3) 
(A20a) 

{Nkgi) + 
(A20b) 



Hence, in order to obtain the unknown expectation values 
(Nkgi), we have to solve the linear system 

{di2k} = -Bi{Nkgi) ^C{Nkg3), (A21a) 
{d23k) - C{Nt,gi)+B3{Nt,g3), (A21b) 



with Ihs. given by Eqs. (jA19|) . This is easily done and 
we may summarize the result for the mean velocity: 

{Rk{t)) = ^ (Nkgi) + ^ {Nkg3), (A22a) 



where 



(Nkgi) = -M^BI0Jd23kl^ (^22b) 

{Nkg3) = (A22C) 

with 

{di2k) = {Nk)di2 + {Nk)di2, (A22d) 

{d23k) = {Nk)d23 + {Nk)d23. (A22e) 

Since the quantities (Nk), (Nk), dap, dap, Aa, Ba and 
C are known we have thus obtained a closed analytical 
result for the mean swimmer velocity within the Oseen 
approximation. Analogous calculations can be performed 
for but do not yield much additional insight (for a 
single swimmer). 

Additional simplifications If the orientation correla- 
tion time Tjs! is larger than the driving period T — 27: /lo 
then 

{di2k) ^ {Nk) rfi2, {d23k) ^ {Nk) d23. (A23) 
In this case, we may simplify 

{Ru{t)) -V{t) {Nk{t)), (A24a) 

where 



V(f\ ^ ( -^3 di2 + Cd23 

^ ' si B1B3 - C2 



A3 / Cdi2+Bi d23 

3 i B,B3-C^ 



(A24b) 



is a periodic function, V{t) = V{t + T^^). Since we as- 
sumed tn ^ T^, we can achieve further simplification 
by replacing V{t) with its stroke-average 



so that 



V := y ds V(s), 



(i?fe(t))^V(7Vfc(t)). 



(A25) 



(A26) 
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For example, when considering equal-sized beads with 
Ua — a and £ 3> max{a, A12, A23}, then 

V[t) = J COS[Ujt + ipi2) H J COS{ujt + ip23) + V, 

where 

V=^ac.(^)sinA^ (A27) 

and Aip :— ipi2 — (p23j and higher order terms have been 
neglected. 



b. Spatial diffusion 

Using the result for dR{t) from above, we may rewrite 
the mean square displacement Dji{t) as 



'Dn{t) ■■= mt)-R{0)r) 
ft ft 
dRk{u) / dRk{s 
/o Jo 

ds / du 
10 Jo 



Ai A3 
-fNkgi + ^Nkg3 



^1 AT ^ ^3 

-J- Nkdi + — Nkg3 

Q X] y([C(afe)(7")]sdB(^„)(s) X 
<y,a' — l 

a'k){'y'n' 

)]„dB(y„,)(u)).(A28) 

Here, we have again used that (/ f{Xa) di3(^„)(i)) = 
holds for Ito integrals. We consider the two remaining 
integrals in Eq. (jA28p separately, starting with the second 
one. We find 



a,a' — l 

a' k)('y'7i')]u 

dB(^,„,)(u)) 

2 . ^ /■* 

a,a' — l ^ 



AAA 

71 72 73 



kB7 ^ (l-^„«0 rd.s?{f„,)(„,,), 

a,a' — l 



where 



arO 

•^(ak){a'k) 



Skk + NkNk 
87r/i daa' 



1 



27r/i daa' 



(A29) 



For spherical particles we have 7^ = 6nfiaa and, there- 
fore. 



1 knJ ( 1 



9 

2kB'y 

9 TT/i 



1 

a2 



t 

V"12 



t + 



1 



1 

d23 



1 

di3 



For £ 3> maxjAct/?}, the integrand in the second line can 
be approximated by 5/{2£) yielding 



^„ , , 1 fesT / 1 1 1 , 



03 



(A30) 



It remains to determine the first (double) integral in 
Eq. ((X28)) . reading 



-DUt) := Tds fdu 
Jo Jo 



-J ^k9i + -J- Nkgs 



Ai 



Nk9i 



A3 



Nkg3 



The subscripts indicate the time arguments in the brack- 
eted expressions, respectively. Upon recalling that 
F{ai) = gaNi is the internal force acting on sphere a, 
we see that the contribution D1^{t) is essentially deter- 
mined by the force-force correlation functions. However, 
instead of calculating these correlation functions exactly, 
we may approximate, for tat ^ T^, the integrand by [cf. 
Eqs. (lIlTl) and (UMll ] 



D^V^ {Nk{s)Nk{u)), 



where V is the stroke-averaged velocity of the correspond- 
ing deterministic swimmer, cf. Eqs. (|A24[) - (|A26[) . Adopt- 
ing this approximation we find 

D%it) ^ f ds f du {Nkis)Nk{u)) 
Jo Jo 

= ds du exp(— |m — s|/rjv) 
Jo Jo 

and thus the final result 

Dfl(t) = DP (t) + V^^it) 

IfcsT/l 1 1 5\ 

- Q— - + - + - + 7 * + 
9 TT/i \ ai a2 03 I J 



2V^TN [t-Hrjv(e-*/^" -1)]. 



(A31) 



The first part represents passive (thermal) diffusion, the 
second part is due to active swimming (note that tn is 
temperature dependent as well). 
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